Current-induced synchronized magnetization reversal of two-body Stoner particles 

with dipolar interaction 
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We investigate magnetization reversal of two-body uniaxial Stoner particles, by injecting spin- 
polarized current through a spin-valve structure. The two-body Stoner particles perform synchro- 
nized dynamics and can act as an information bit in computer technology. In the presence of 
magnetic dipole-dipole interaction (DDI) between the two particles, the critical switching current I c 
for reversing the two dipoles is analytically obtained and numerically verified in two typical geomet- 
ric configurations. I c bifurcates at a critical DDI strength, where I c can be decreased to about 70% 
of the usual value without DDI. Moreover, we also numerically investigate the magnetic hysteresis 
loop, magnetization self-precession, reversal time and the synchronization stability phase diagram 
for the two-body system in the synchronized dynamics regime. 

PACS numbers: 75.60.Jk, 75.75.-c, 85.75.-d 



I. INTRODUCTION 

Recent advances on nanomagnetism technology allow 
to fabricate Stoner particles (single-domain magnetic 
nanoparticles due to strong exchange interaction)^— , 
which has attracted much attention both from a funda- 
mental point of view and potential applications in infor- 
mation industry^ 6 -. Magnetization dynamics of a single 
Stoner particle has been extensively studied, either by 
static or time-dependent magnetic fields^—, or by spin- 
polarized current approaches^—. The current-induced 
magnetization switching attracts much interest due to 
the locality and convenient controllability of currents. 
However, the high critical switching current I c (or den- 
sity) hinders practical applications^. Many efforts were 
made for lowering the I Cl for instance, designing an op- 
timized spin current pattern^, using pure spin currents! 
and thermal activatio n 25 ' 26 . Recently lower I c has also 
been proposed in ferromagnetic semiconductors^^. 

Moreover, since magnetic nanoparticles are actually 
fabricated in arrays^—, the dipole-dipole interaction 
(DDI) between them will be important for the magnetic 
switching behavior. Indeed there were several theoretical 
studies on systems of two Stoner particles, the simplest 
case to investigate the DDI effect, concluding the dipolar 
interaction can assist the magnetic switching22r— . Re- 
cently a collaboration including the present authors re- 
visited the two-body Stoner particles problem and pro- 
posed a novel technological perspective by showing that, 
in a synchronized dynamics mode (i.e. both the magnetic 
vectors of the particles runs synchronously) , the critical 
switching field in presence of magnetic DDI can be dra- 
matically lowered. This implies even zero-field switch- 
ing to be achieved by appropriately engineering the DDI 
strength by adjusting the interparticle distance—. In the 
present paper, we extend these studies by investigating 
current-induced synchronized magnetization switching of 
such the two-body Stoner particles system on injecting 
a spin-polarized current through a spin-valve-like struc- 
ture. As a result, the critical switching current I c (or 



its density) can be lowered to about 70% of the usual 
value without DDI around a critical DDI strength. This 
DDI strength corresponds to a distinguished interparti- 
cle distance and has the same value as in the previously 
studied case of field-induced switching^ 6 -. Furthermore, 
I c bifurcates at the critical DDI strength by a square- 
root-like behavior on the DDI parameter. Moreover, we 
have also numerically investigated issues including the 
magnetic hysteresis loop, magnetization self-precession, 
reversal time, and the synchronization stability phase di- 
agram for the two-body system in presence of the DDI. 
These investigations revealed many novel promising phe- 
nomena regarding future device applications. This paper 
is organized as follows: In Sec. II the two typical geo- 
metrical configurations of the two-body system in a spin- 
valve-like setup are introduced along with the equations 
governing current-induced magnetization dynamics. Our 
analytical and numerical results are presented in Sec. Ill 
and IV, respectively. We close with a discussion and con- 
cluding remarks in Sec. V. 



II. MODEL 

We consider the spin-valve structure schematically 
shown in Fig. [TJ For simplicity, we consider two par- 
ticles are in the same sphere forms. Two typical ge- 
ometric configurations of the two-body system are in- 
vestigated, where the unit vector n along the line con- 
necting the two particles is either perpendicular or par- 
allel to the magnetic easy axes (EA), assumed along 
z-axis, referred to as PERP configuration and PARA 
configuration, respectively. That is, if we define h = 
(sin 9 n cos <ftn, sin^n sin </>„, cos#„) in the usual spherical 
coordinates where z-axis is the pole axis, thus 9 n = n/2 
for PERP configuration and 9 n — 0, tt for PARA (the 
two values are equivalent since DDI is invariant under 
h M> — n). Moreover, without loss of generality, one can 
let <fi n = for both configurations, i.e., n = x in PERP 
and h = z in PARA (see Fig. p. 
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FIG. 1: A schematic diagram of two typical geometric con- 
figurations (PERP and PARA) in spin-valve structures and 
the corresponding coordinate axis orientations. Polarizer is 
a thick pinned ferromagnetic layer where s denotes its polar- 
ized direction and NM denotes nonmagnetic metallic layer. 
EA denotes the easy axis and e denotes the electrons flow. 



The magnetization dynamics of two-body Stoner par- 
ticles under a spin-polarized current is governed by the 
modified Landau-Lifshitz-Gilbert (LLG) equation 37 with 
the spin-transfer torque (STT) ter m 19 i 20 , 



X hi + arhi x rhi + T sl 



(1) 



Here rhi = Mi/M s is the normalized magnetization vec- 
tor of the ith-particle (i = 1,2), M s = |Mj| is the mag- 
netization saturation of either particle, and a is the phe- 
nomenological Gilbert damping coefficient. We are as- 
suming the two particles are completely identical with all 
the same concerned parameters for simplicity. The unit 
of time is set to be (|7|Af s ) _1 where 7 is the gyromagnetic 
ratio. The total effective field h\ on each particle comes 
from the variational derivative of the total magnetic en- 
ergy with respect to magnetization, h\ = —5E/5r%i. For 
the concerned system, the energy E per particle volume 
V (in units of ^lqM^; /_i : vacuum permeability) can be 
expressed as, 



E=- 



i=l,2 



kmf z + rj[ihi ■ 777,2 — 3(mi • h){fn,2 • n)]. (2) 



Here the uniaxial magnetic anisotropic parameter k sum- 
marizes both shape and magnetocrystalline contributions 
to the magnetic anisotropy along the easy axis and we 
only consider no external magnetic field case. The pa- 
rameter rj = is a geometric factor characterizing 
the DDI strength, where d is the fixed distance between 
the two particles. We also omit the exchange inter- 
action energy between two particles as the same rea- 
sons discussed in Ref. l3g. The LLG equation can be 
written in the usual spherical coordinates, by putting 

rhi = ( sm @i cos 4>i > sin @i sm 4>i j cos @i ) ■ 

The important STT term T sft = —ajrhi x (m, x s), 
(i = 1,2), taking the Slonczewski's formic where ai — 
2noeM' 1 v nas a (normalized) magnetic field dimension. 
We omit the field-like STT terrn^ because it is usually 
smaller than the Slonczewski's term. Furthermore, it has 



the magnetic field character whose role for the two-body 
system was already discussed previously- 3 ^. Here (among 
standard notation) I and P are the magnitude and spin 
polarization degree of the current, respectively, s is the 
unit vector along the current polarization direction (see 
Fig. [T]) which is characterized by two angles 9 S and <p s 
in the spherical coordinates. Since a/ is proportional to 
the current I (assuming other parameters are constant), 
we will refer ai as to the switching current in the follow- 
ing description. For instance, a/ jC stands for the critical 
switching current. 



III. ANALYTICAL RESULTS 

In the following we will concentrate on the synchro- 
nized motion mode of the two dipoles in the current- 
induced switching. That is to say, both magnetiza- 
tion vectors remain in parallel throughout the motion, 
Q\ = 62 = 9, and (pi = 4>2 = <p- Analytical results can 
be obtained in such motion and the stability of the syn- 
chronization will be numerically verified later. Thus, the 
two-body Stoner particles promise to play the role of an 
information bit in computer technology. The modified 
LLG equations ([T]) in spherical coordinates are a system 
of coupled nonlinear equations and read 



V + a sm ( 
a9 — sin ( 



COS tpg 

-ai — — — - 377 cos ip n sin 9 n sin (p, 
o9 

-ai sin 9 S sm(<p — <p s ) — k sin 29 
3rj d cos 2 ip n 
+ Y 09 ' 



(3) 



Here ip n /s is the angle between rh and fi/s such that 
cosi/vi/s = cos#cos# n / s + sin 6> sin (9„/ s cos(0 — <p n / s ), 
(Note that we define <p n — 0). Eq. © is the starting 
point for the numerical simulations. 

In order to find the critical switching current a/ jC , one 
should find the fixed points (FPs) of the nonlinear equa- 
tions ([3]) and then analyze their stability conditions 38 . 
Therefore, we first examine the FP condition for the two 
poles 9 = 0, 7T. From Eq. ©, we obtain 

(3/2)77 sin 26 n (a cos <p qp sin (p) 

+ aj sin 6 S [a sm((p — <p s ) ± cos(0 — (p s )] = 0, (4) 

where =p or ± corresponds to 6 — and 6 — it, respec- 
tively. Due to the uncertainty of (p at poles, we conclude 
only sin2(9„ = and sm# s = such that the two poles 
become the possible FPs. Thus, we only need to focus 
on the simple cases of 6 S = 0,ir and 9 n = Q,tt,tt/2 (i.e. 
PARA and PERP configurations). In the following, we 
let 9 S = 7T, which is the equivalent situation for 9 S = 
by inverting the current flow direction. 

Now let us first recover the usual case without DDI 
(i.e. rj = 0), where PERP and PARA configurations are 
indistinguishable. The motion of 9 is decoupled to <p, 



T9 = a j sin 9 — ak sin 29, 



(5) 
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where we define T = 1 + a 2 . The critical switching cur- 
rent a/ ;C can be obtained by analyzing the linear stability 
conditions at poles 3 -^. After some algebra one can easily 
find when aj > 2ak the FP of 9 = becomes unstable, 
but 9 = 7r is stable, which implies that the magnetization 
reversal switches on and a/ jC reads 

a° Ic = 2ak, (6) 

where the superscript index denotes the zero DDI case. 
We also like to point out that 2k is the usual Stoner- 
Wohlfarth limiti in the field reversal case, namely the 
barrier height between the two minima at poles in the 
energy landscape. 

We now turn to the case with DDI (r; ^ 0). Let us first 
write down the magnetostatic energy (without external 
field) under the synchronized motion mode: 

E = -2fc cos 2 6> + 7/(1 - 3cos 2 i/>„), (7) 

which immediately leads to the equilibrium (minimal) 
energy for the PARA and PERP configurations as 

e para =v _(2k + 3rj) cos 2 9, 

E PERP = _ 2l] + (3?? _ 2fc j CQS 2 q (g) 

Fig. [2] shows that the essential differences in the en- 
ergy landscapes for the PARA and PERP configurations. 
Figs. Ufa)- (c) show the energy variation in PARA config- 
uration with the increase of the DDI strength 77. The two 
edges of 9 — 0, n should be understood as the two poles 
in a sphere. In PARA cases, the energy diagrams have 
no essential geometrical shape change with variation of 77. 
The two poles (9 = 0, tt) are always the energy minima, 
and the in-between barrier (at 9 = 7r/2) height increases 
with 77, expecting a higher critical switching current (or 
field) in the PARA case. 

However, in the PERP configurations as shown in 
Figs. Hld)-(f), the energy landscapes exhibit an essen- 
tial geometrical change with the increase of 77. Re- 
markably, from Eq. (|5J|, a critical DDI strength r\ c = 
2k J '3 exists, which is the same value both in the cur- 
rent and field switching cases. Furthermore, r\ c is 
equivalent for the interparticle distance being d c — 
{S^oM^V/SttK) 1 ^, where K = k(j, M% is the standard 
anisotropy coefficient 3 ^. At the critical distance, the en- 
ergy diagram has two flat paths along <j> = 0, <j> [see 
Fig. H[e)], which means any position of 9 can be equi- 
libriums. While on the both sides of the critical t] c value, 
the system synchronized ground states (i.e. energy min- 
imum) change. When the interparticle distance d > d c 
(t) < Vc), the ground states are at two poles of 9 = 0, tt 
[Fig. [H(d)] while when d < d c (77 > r] c ), the ground states 
are at another two positions, 9 = n/2 and <f> — 0, tt, along 
the hard x-axis. 

For the PARA configuration, one can find the decou- 
pled equation of 9 to </>, 

= ai sin 9-a(k + 3r)/2) sin20. (9) 



(a) PARA(I7=0.1) (d) PERP<^=0.1) 




(b) PARA(^=l/3) (e) PERP(j;=l/3) 
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(c) PARA(q=0.9) (f ) PERP(i;=0.9) 




FIG. 2: (color online.) The three-dimensional energy di- 
agrams for the PARA and PERP configurations, (a)-(c): 
PARA with 77 = 0.1,1/3,0.9; (d)-(f): PERP with 77 = 
0.1, 1/3, 0.9. The system parameter k = 0.5. 

With the same linear stability analysis as before, one can 
derive that when a/ > a(2k + 3rj), the FP of 9 = 
becomes unstable, but 9 = tt is still stable. Hence, the 
critical switching current in the PARA case reads, 

alf RA = a(2k + 377). (10) 

Thus, af ARA is always higher than aPj for 77 ^ 0, as 
is expected that the barrier height between two energy 
minima is always increased with increasing 77. 

The PERP configuration is much more interesting and 
exhibits novel physical behavior due to the geometrical 
structure change in the energy diagrams as shown in 
Figs.[UJd)-(f), in which a novel zero-field switching mech- 
anism was proposed 3 ^. In PERP, the nonlinear Eq. ^ 
cannot be decoupled any more for 9 and 0, 

T(9 = ai sin 9 — ak sin 29 

+ (3/2)7i(asin26>cos 2 <j> - sin sin 20), (11) 
T(f> = aaj + 2k cos 9 — 3r/ cos 9 cos 2 (f> 

- (3/2)a77sin20. (12) 

Again according to the linear stability analysis approach, 
the linearized matrices at the FPs of 9 — 0,tt can be 
found as A\g =0 ^ = diag(Af 1; Af 2 ) with 

Af-L = ±a/ -2afc + 377(0; cos 2 (sin 20)/2) (13) 
A± 2 = 377(±sin2(/5-acos20) (14) 
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where diag denotes the diagonal matrix and the upper 
(lower) sign refers to 9 = (9 = ir). After some algebra, 
one finds that for a/ > a(2k — 3r]/2), 9 = becomes 
unstable but 9 = ir is still stable, such that the critical 
switching current reads, 

a PERP,l = a{2k _ 3rj/2)) {71<lh) (15) 
PERP 1 

The obtained a,j ' can only be applicable for the 
regime of r\ < r\ c because the two-body system in the 
PERP case changes its energy equilibriums for rj > r) c , 
where 9 = n/2 may become the FP, as shown in energy 
diagram Figs.[Hd)-(f). Thus, one requires to analyze the 
stability condition around the new FP for the 7/ > r/ c 
regime, which can be obtained under the current case as: 

9* = tt/2, 4>* = (sin~ 1 (2a / /3r/)) /2. (16) 

Then, the linearized matrix at this FP is 

_ (2ok - Za-q cos 2 </>* -3?7Cos2</>* \ , . 
A\e\r - ^ 3 7?cos 2^*_ 2 fe -3ar) cos2(/>* J ■ [U) 

After some algebra, we can find that when 4a 2 - + (4fc — 
3r]) 2 > 9r] 2 , the FP becomes unstable, but the pole 9 = tt 
is still stable. Thus the critical switching current within 
the r/ > r\ c regime reads 

a P f Pp2 = y/2k(3r,-2k), (77 > Vc ). (18) 

So far we obtained the analytical solutions of the criti- 
cal switching current in both the PARA and PERP con- 
figurations, aff RA , af^ RP ' 1 and a PPPp2 , which are 
plotted by the solid, dashed and dot-dashed lines in 
Fig. 131 respectively. The system parameters k = 0.5 
and the damping a = 0.1. In principle, the minimum 
switching current in PERP should occur at the critical 
DDI value 7? c , 

aZ n = ak = a° IiC /2, (19) 

which is obtained by using Eq. (|15j) at rj c and is a half of 
the value without DDI. However, as we will show later 
[see Fig. |SJa)], there is a very small region around r\ c 
above the crossed two solutions Eq. (fTS")) and Eq. (Ti"8l) . 
where the magnetization dynamics is not static but shows 
self-precession. Thus, in practice, we numerically found 
that the minimal critical switching current a™" for a 
static magnetization reversal can be lowered as about 
70% of the usual value a® c without DDI, which will be 
more discussed later. 



IV. NUMERICAL RESULTS 

In order to complement and quantitatively support our 
previous discussions, we will now numerically verify our 
analytical results and investigate the reversal phenomena 
of two-body Stoner particles system in presence of mag- 
netic dipolar interaction under a spin-polarized current. 
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FIG. 3: (color online.) The (normalized) critical switching 
current a/ iC versus the DDI strength r\ for the PERP and 
PARA configurations. The lines show analytical results and 
the circles and squares show numerical results. The system 
parameters are k = 0.5 and a = 0.1. 

We solved the Eqs. (|3|) by implementing a fourth-order 
Runge-Kutta method. We consider the range of the DDI 
strength to be < rj < 1. In this case, rj = corresponds 
to the limiting case of two particles infinitely apart from 
one another i.e. d — > oo. Large i] may be realized by fab- 
ricating magnetic nanoparticles of ellipsoidal shapes al- 
lowing for a closer proximity. Throughout the numerical 
calculations, we use the damping parameter of a = 0.1. 
In Fig. [3j we first compare the numerical simulations 
of the critical switching current a/ jC in both PERP and 
PARA configurations (circles and squares) and show an 
excellent agreement to our analytical solutions: Eq. (|10l) . 
Eq. ([15}) and Eq. (fTS)) . The critical switching current 
in PARA is always higher than the value without DDI 
(at rj = 0), while in PERP it bifurcates at the critical 
DDI value (at 77 = 1/3) with a square-root dependence 
behavior. 



A. Hysteresis loop 

We then, with the numerical simulations, study the 
magnetic hysteresis loops under the current-induced syn- 
chronized switching. Fig. Hta) first shows the hysteresis 
loops at different DDI strengthes for the PARA config- 
uration. The increase of rj will enlarge the loop sizes. 
Fig. HJb) shows the hysteresis loops for the PERP con- 
figuration for r] < r) c . The increase of rj will reduce the 
loop sizes, in contrast with the PARA case. In particular, 
when ry approaches closely to r) c , magnetic self-precession 
phenomenon in current switching can be observed, which 
will be discussed in detail in next subsection. For exam- 
ple, there are very narrow regimes at rj = 0.31, which is 
enlarged in the insets of panel (b), showing the m z oscilla- 
tion amplitudes. Fig. |4jc) shows the hysteresis behavior 
for the PERP configuration for 77 > r\ c . It is interesting 
to note that there is no any loops but a reversing tri-state 
occurring. At last, as a comparison, we plot panel (d) to 
show the hysteresis behaviors in the field switching case 
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for the PERP configuration, where h is the normalized 
external field by the unit of magnetization saturation M s . 
When rj < r) c a clear hysteresis loop occurs while when 
f] > Vc there is neither hysteresis loops but the magne- 
tization m z is linearly stable with h, different from the 
tri-state phenomenon in the current driven case. 



1.0- 
0.5 

E o.o- 

-0.5 

-1.0- 



- ii=0.2 
-,=0. 

- 

-0.4 -0.2 0.0 0.2 0.4 0.6 

a, r 




-0.07 -0.06 -0.05 — 

-0.4 -o.2 o!o 0^2 04 




PERP (field) 

11=0.1 

— o ii=0.5 



! Hi 



(d)- 



-1 

h[MJ 



FIG. 4: (color online.) Hysteresis loops, (a) PARA config- 
uration, (b) PERP configuration when r\ < r\ c — 1/3. The 
insets show the enlargement of the magnetic self-precession 
regimes, where m z is proportional to the oscillation ampli- 
tude, (c) PERP configuration when r) > r/ c = 1/3. (d) PERP 
configuration in the field switching case. Other parameters 
are the same as Fig. [3] 



B. Self-precession 

With the numerical simulations, we observed in the 
PERP configuration, there exists a small region above 
the crossed two solutions Eq. (fl"5j) and Eq. (TT5|) around 
the critical DDI value r\ c . The magnetization vectors in 
this region perform the STT induced self-precession (or 
self-oscillation) phenomenon—-.^. In Fig. EJa), we plot 
the self-precession phase diagram in the parameter rj — aj 
space. The plotted color is proportional to the magnetic 
self-oscillation amplitude for m z (t). The self-precession 
occupies only a very small partition in the parameter 
space for the current-induced synchronized dynamics of 
two-body systems, which is consistent with the previous 
studies for single biaxial magnetic nanoparticle without 
DDI 38 . Fig.[5]Jb) shows an example of the self-oscillation 
of m z (t) in time domain at 77 = 0.31 and aj = 0.06. The 
inset shows the spacial trajectories of the magnetization 
vector and its projections in m x — m y — m z coordinates. 
We also like to point out that the practical critical switch- 
ing current a™™ at r\ = i] c cannot be lowered to a half of 
that without DDI due to the self-precession phenomenon. 
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FIG. 5: (color online.) (a) The phase diagram in the PERP 
configuration of the magnetization self-precession in the pa- 
rameter r\ — ai plane, where the color is proportional to the 
self-oscillation amplitude and the two lines show the analyt- 
ical solutions of the critical switching current, (b) The self- 
oscillation curve of m z (t) at rj — 0.31 and ai — 0.06. The 
inset shows the spacial trajectories of the magnetization vec- 
tor and its projections. Other parameters are the same as 
Fig. O 



C. Reversal time 

The numerical results for the reversal time T r versus 
the switching spin-polarized current magnitude aj at dif- 
ferent DDI strength for both PARA and PERP configu- 
rations are shown in Fig.|5[a) and (b). The reversal time 
is numerically determined from m z = +1 to m z — — 1 for 
either particle, which depends on the initial directions of 
the magnetization vectors due to the metastability of the 
polesis. Thus, in practice, we deviate the initial angles 
6*i (0) = 62(0) — 0.001 in order to find a finite reversal 
time. The same initial values for either particle is of 
importance for the stability of the synchronized motion 
mode, as we will discussed in next subsection. Fig. [BJa) 
and (b) show that, in both the PARA and PERP cases, 
the reversal time curves with current have similar trends, 
approximately inversely proportional to the current mag- 
nitude a 1 for large currents. In the PERP configuration, 
the T r behavior is different from that in the field switch- 
ing case when the DDI parameter approaches its critical 
value. A substantially shorter reversal time can be found 
around zero field regime in the field case 3 -^ when rj m r\ c . 
In contrast, the reversal time curve is non-smooth in the 
current case when rj is around its critical value, as shown 
in the inset of Fig. [BJb). 



D. Stability phase diagram 

As an important issue, we now examine the stability 
of the synchronized magnetization dynamics scenario in 
the PERP configuration. We examine the stability from 
two aspects. One is that the difference in initial angles of 
the two magnetization vectors may result in unstablity 
for the synchronization. The other is the tiny difference 
in two particles' material parameters may spoil the syn- 
chronization. 

Let us first numerically examine the average stable 
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FIG. 6: (color online.) Reversal time T r versus the cur- 
rent magnitude ai. (a) PARA configuration: different DDI 
strength r\ — 0,0.1,0.5,1; (b) PERP configuration: r\ — 
0,0.1,0.2,0.33,0.4,0.7. The inset is an enlargement for 77 = 
0.33 where non-smooth T r is observed. The connected lines 
are guides to eyes. Other parameters are the same as Fig. [3] 



magnetization (m\ z + m2 Z )/2 (the value —1 means stable 
while unstable) with the dependence of the deviation 
angle S between the initial directions of the two magne- 
tization vectors. In detail, we fixed the initial direction 
of one particle to be along the z-axis i.e. 61 = 0.001 and 
changed the other one as 82 = 0.001 + 6 where 0.001 is 
for eliminating the singularity in poles. We can plot the 
stability phase diagram in the ai — r\ space, as shown 
in Fig. [Tfa) for the current case, and plot Fig. [Tfb) for 
the field switching case as a comparison in the h — r) 
space. The color difference is proportional to maximal 
8 m in unit of degrees for guaranteeing the stability where 
the gray color denotes S m > 1° in (a) and S m > 10° in 
(b). That is to say, when 8 > 5 m , the final average mag- 
netization is found to be 0, implying the synchronized 
dynamics is destroyed. One can observe two points from 
Fig. [TJa) and (b): (1) In both current and field switch- 
ing cases, the important stability regimes are located just 
above the critical switching current or field lines shown 
in dashed or dot-dashed lines in corresponding figures, 
especially for 77 < rj c . "Islands" appear for such regions. 
(2) It has wider initial deviation angles limit 5 m for the 
synchronization stability in the field case (< 10°) than in 
the current case (< 1°). We also plot an example shown 
in Fig. [8] in the case of 77 = 0.23 and ai = 0.07, which 
is 70% of the usual switching current =0.1 without 
DDI. Fig. EJa) shows the time evolution of all the com- 
ponents of the magnetization vector. Fig.^b) shows the 
maximal initial deviation angle S rn » 1.2° for the syn- 
chronization stability. Consequently, we remark that the 
practical critical switching current a™" for a static mag- 
netization reversal can be lowered as about 70% of the 
usual value a® c . 

On the other hand, for two particles with slightly dif- 
ferent, for instance, anisotropy Ki, volume Vi, satura- 
tion magnetization M Syi (i — 1,2), as discussed before^, 
the critical value of DDI becomes r\ c = (ki + fc2)/3 = 

3^m^mmv^+v 2 ) and the effective neld experienced by 
each particle due to the external field or the STT here will 
be slightly different. Hence, we also numerically examine 
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FIG. 7: (color online.) Synchronization stability phase dia- 
grams in the PERP configuration for the initial direction de- 
viation 5 of the magnetization vectors in current case (a); in 
field case (b) . The color difference is proportional to maximal 
stability 5 m in unit of degrees where the gray color denotes 
S > 1° in (a) and S > 10° in (b). Synchronization stability 
phase diagrams for different experienced Sai in current case 
(c); Sh in field case (d). The color difference is proportional to 
the variation width Saf or 5h w where the gray color denotes 
Saf > 0.2 and Sh w > 0.02. Other parameters are the same 
as Fig. El 
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FIG. 8: (a) The magnetization vector components versus time 
at r\ — 0.23 and ai — 0.07. (b) The average stable magnetiza- 
tion (miz +m2z)/2 versus the initial angle deviation 5, where 
S m ~ 1.2° is observed, (c) (mi z + wi2z)/2 versus the current 
deviation Sai, where 8a™ m 0.04 is found. Other parameters 
are the same as Fig. [3] 
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the synchronization stability due to such difference. In 
detail we fixed the current torque on one particle as a\j 
and changed the other one as 0,2, i — 01,1 + Saj where 
8a 1 denotes a small deviation for different currents. The 
field case is the same except introducing a small devia- 
tion 5h = I12 — hi for different external fields. We thus 
plot the stability phase diagrams in the aj — 77 space, as 
shown in Fig. [Tfc) for the current case, and Fig. EJd) for 
the field case as a comparison in the h — rj space. The 
color difference is now proportional to the variation width 
of 8af and Sh w , in which the final average magnetization 
is found to be —1 for stability. The gray color denotes 
Saf > 0.2(M S ) in (c) and 5h w > 0.02(M S ) in (d). From 
the figures, one can observe the stability "islands" struc- 
tures above the critical switching current or field lines 
become smeary now, and the current switching case has 
a wider stability window due to the parameter differences 
of the two particles than that in the field case. Fig. \%lc) 
shows an example for current difference on two particles 
at the case of i] = 0.23 and a/ = 0.07, where the width 
Saf * 0.04 for the average stable magnetization equaling 
— 1 to sustain the synchronized motion mode. 

V. DISCUSSION AND CONCLUSION 

We first like to compare our results with a concrete 
magnetic material such as cobalt (Co) particles. The 
standard data is M s — 1400kA/m, uniaxial strength K = 
10 5 J/m 3 , and a = 0.1 13 . Thus k = K/(fi M^) = 0.04 
such that the critical DDI strength r\ c — 0.027. If we con- 
sider two spherical particles with the radius r so that the 
DDI parameter 77 = r 3 /(3d 3 ). Thus the critical DDI is 
reached at d c — 2.3r. Assuming r = lOnm and spin po- 
larization P = 0.4, the critical switching "effective field" 
without DDI is a° I c = (2ak)M s = 140Oe, which trans- 
lates the critical switching current density to be about 
jo = 3 x 10 8 A/cm 2 . In our model with DDI we argue that 
the critical switching current can be lowed to be 70% of 



the original value, i.e. about 2 x 10 8 A/cm 2 when the two 
particles are engineered to be located near the critical 
distance d c = 23nm. Also in the case of Co, the unit of 
time is (|7|Af s ) _1 = 3.23ps. From Fig. [f^b), the rever- 
sal time is infinitely long at the critical switching current 
point. Increasing the switching current will drop the re- 
versal time inversely. More importantly, the key issue 
for the technological aim in our two-body Stoner parti- 
cles system is to maintain the synchronized motion mode 
against large deviations from such as initial conditions 
and/or material parameters. The stability region from 
the initial angle deviations in the current case (about 1°) 
is much smaller than that in the field case (about 10°), 
thus how to enhance the synchronization stability might 
be most interesting for future studies. Moreover whether 
there exists a zero switching current in contrast with the 
zero-field switching case with the aid of DDI is also an 
interesting issue. 

In conclusion, we have investigated the magnetization 
reversal of two-body uniaxial Stoner particles in a stable 
synchronized motion mode, by injecting a spin-polarized 
current through a spin- valve like structure. In presence of 
magnetic dipolar interaction, the critical switching cur- 
rent for reversing the two dipoles is analytically obtained 
and numerically verified in two typical geometric PERP 
and PARA configurations. In the interesting PERP con- 
figuration, the critical switching current bifurcates at a 
critical DDI strength with a square-root behavior, near 
where it can be lowered to about 70% of the usual value 
without DDI. Moreover, we also numerically investigate 
the current-induced magnetization hysteresis loops, mag- 
netic self-precession phenomenon, reversal time and the 
synchronization stability phase diagram in the two-body 
system, which shows interesting predictions and is ex- 
pected to be useful for future device applications. 

Z.Z.S. thanks the Alexander von Humboldt Founda- 
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ported by Deutsche Forschungsgemcinschaft via SFB 
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